

NASA TECHNICAL NOTE 





I 





,^(NASA-TN-b-7410) A RAPID NDMERICAL N74-16703^ 

SOIOTIOH TO SUBSONIC FLOW OVER PLANAR 
;AND AXISYMMETRIC PROFILES AT AN ANGLE OF 

.ATTACK OF 0 DEG (NASA) p HC $3.25 Unclas 

CSCL OlA HI/01 30371 J 



A RAPID NUMERICAL SOLUTION 
TO SUBSONIC FLOW OVER PLANAR 
AND AXISYMMETRIC PROFILES 
AT AN ANGLE OF ATTACK OF 0“ 


by Vincent R. Mascitti 

Langley Research Center 
Hampton, Va. 23665 


NATIONAL AERONAUTICS ANO SPACE ADMINISTRATION • WASHINGTON, D. C. • MARCH 1974 




For sale by the National Technical Information Service, Springfield, Virginia 22151 


1. Report No. 


2. Government Accession No. 


3. 


Recipient's Catalog No. 


NASA TND-7410 


4, 


Title and Subtitle 

A RAPID NUMERICAL SOLUTION TO SUBSONIC FLOW OVER 
PLANAR AND AXISYMMETRIC PROFILES AT AN ANGLE OF 
ATTACK OF 0° 


5. Report Date 

March 1974 

6. Performing Organization Code 


7. Authof^s) 


8. Performing Organization 


Report No. 


Vincent R. Mascitti 


L-8904 

10. Work Unit No. 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


' ' 501-24-05-01 

11. Contractor Grant No. 


13. Type of Report and Period Covered 

Technical Note 

14. Sponsoring Agency Code 


16. Abstract 

A rapid numerical solution is presented in this paper for the incompressible flow over 
thin planar and axisym metric profiles at an angle of attack of 0^. The method uses a finite - 
difference field solution to the governing equation with a Gauss-Seidel successive overrelaxa- 
tion scheme. However, the use of a simple Cartesian grid system restricts this method to 
slender profiles. Results are presented for a cambered airfoil, airfoil in wall effect (two- 
dimensional flowthrough inlet), body of revolution, and flowthrough nacelle. 

A computer program is presented which can be used for any of the previously mentioned 
cases with simple input changes. Results for compressible flow are available with the"use of 
the appropriate two-dimensional or axisymmetric compressibility corrections. Computational 
time for a typical field calculation of 3000 grid points and 200 cycles through the field is less 
than 1 minute with less than 50 000 octal storage on the Control Data Corporation 6600 com- 
puting system at the Langley Research Center. 


17. Key Words (Suggested by Author (s)| 

Subsonic flow 
Two-dimensional profiles 
Finite difference 


18. Distribution Statemer^t 

Unclassified - Unlimited 


STAR Category 01 


21. No. of Pages 


22. Price 3.^^ 

■mm 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 









r 


A RAPID NUMERICAL SOLUTION TO SUBSONIC FLOW OVER PLANAR 
AND AXISYMMETRIC PROFILES AT AN ANGLE OF ATTACK OF 0° 

By Vincent R. Mascitti 
Langley Research Center 

SUMMARY 

A rapid numerical solution is presented in this paper for the incompressible flow 
over thin planar and axisymmetric profiles at an angle of attack of 0®. The method uses 
a finite -difference field solution to the governing equation with a Gauss-Seidel successive 
overrelaxation scheme. However, the use of a simple Cartesian grid system restricts 
this method to slender profiles. Results are presented for a cambered airfoil, airfoil in 
wall effect {two-dimensional flowthrough inlet), body of revolution, and flowthrough 
nacelle. 

A computer program is presented which can be used for any of the previously men- 
tioned cases with simple input chajiges. Results for compressible flow are available with 
the use of the appropriate two-dimensional or axisymmetric compressibility corrections. 
Computational time for a typical field calculation of 3000 grid points and 200 cycles 
through the field is less than 1 minute with less than 50 000 octal storage on the Control 
Data Corporation 6600 computing system at the Langley Research Center. 

INTRODUCTION 

In the past decade the advent of the high-speed digital computer and the development 
of numerical techniques such as finite differences has enabled the aeronautical researcher 
to attack many problems which were previously not tractable. Especially in low-speed 
aerodynamics, a diversity of problems which consider nonlinear effects with inviscid 
formulations have been approached. 

Solutions to the subsonic and transonic flow over airfoils are presented in refer- 
ences 1 and 2. The method of reference 1 is formulated in terms of the velocity potential 
and rapid computational times are realized even for the transonic flow case. However, 
extensions of the method to rotational or viscous flow is prohibited. The method of ref- 
erence 2 was formulated in terms of the stream function, but it was made clear that the 
stream function formulation is probably not feasible for transonic flow because of the 
difficulty in evaluating the density function, which becomes double valued in the region of 



ti*3.nsonic flow. Tho tF3.nsonic flow 9.l30ut slondOF bodies of revolution ws.s solved in ref~ 
erence 3 with very good accuracy. Once again the velocity potential formulation was 
utilized. The solution to the transonic nacelle problem was performed in reference 4 by 
using the streamtube curvature approach. 

With these methods available for subsonic and transonic flows, it is appropriate to 
consider a general solution procedure for subsonic flows which encompasses these pro- 
files and is extendable to rotational and viscous flows. 

The purpose of this paper is to present a simple and rapid approach to the general- 
ized two-dimensional (planar or axisymmetric) flow problem requiring the solution of 
elliptic differential equations. The computer program can be used to reproduce classical 
incompressible flow solutions, such as the airfoil and body of revolution. But more 
important, the technique can be used as a springboard for more complex problems in 
two-dimensional subsonic flow. The program is simple enough to insure that the built-in 
algorithm for the present governing equation can be replaced with others corresponding to 
different equations. For example, the rotational flow problem with nonuniform upstream 
velocity can be solved as in reference 5. A simplified Cartesian grid system which 
restricts the calculations to slender sections at an angle of attack of 0° is used. Grid 
refinement in the axial direction is employed to provide a better definition of rapidly 
changing flow properties. A detailed listing of the computer program with operating 
instructions is presented. 

Illustrations of program applications are presented for cambered airfoil and airfoil 
with thickness in wall or ground effect, a parabolic body of revolution, and a camber-line 
nacelle. Results are compared with classical conformal mapping or finite- element sin- 
gularity solutions. A comparison with experimental data for a NACA series nacelle hav- 
ing thickness and camber, at a Mach number of 0.6 and a mass-flow ratio of 0.787 is also 
presented to demonstrate the versatility of the method. 

SYMBOLS 


Cp pressure coefficient 

c chord 

LJ index for axial and lateral field points (see fig. 1) 

K two-dimensional flow index. If = 0, planar; if = 1, axisymmetric 

L matrix size 
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m^niQ mass-flow ratio (the mass flow captured by the profile compared with the 
flow captured in free stream) 

R iteration number 

Uq free -stream velocity 

u,v velocity components in the z- and r-directions, respectively, nondimension- 

alized to free-stream velocity 

2,r coordinates (see fig. 1) 

6 angle corresponding to average trailing-ec^e slope (positive clockwise, see 

fig. 1) 

V overrelaxation factor 

]// stream function 

upstream value of stream function at profile radius 
METHOD OF SOLUTION 

The governing equation for steady, irrotational, incompressible flow in two dimen- 
sions is given in reference 6 as 

+ = o (1) 

ar2 r 8r 3^2 

If K = 0, the flow is planar. If K = 1, the flow is axisymmetric. The velocity compo- 
nents are given by: 



( 2 ) 
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In the present finite -difference field solution, it was assumed that an infinite field can be 
represented by a finite field provided that it is sufficiently large compared with the body. 
Furthermore, along the boundaries of the finite field, the boundary condition at infinity 
was applied; that is, v = 0 and u = 1. Therefore 

cc i-K+1 

^boundary 

The boundary condition on the surface of the profile is satisfied by 
'^surface “ Constant 


The present method accounts for unequal spacing around a grid point, but for clarity 
the following equations and algorithms are presented for equal spacing. The finite - 
difference algorithms for unequal spacing can be readily derived and are presented in 
appendix A. 

The second-order finite- difference algorithm for equations (1) and (2) with equal 
grid spacing is 


^(■'' 1 , 1.1 - * 'I'i.i-i 


K 

2rj 







(3) 



( 4 ) 


For a given number of grid points L, equation (3) yields L simultaneous linear 
equations which can be solved by the inversion of an L x L matrix. However, the result- 
ing influence matrix is strongly diagonal and suggests a solution by iteration. 

A possible method is the Gauss- Seidel iteration scheme given by 
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With this method the calculation marches (fig. 1) from lower left to upper right and 
continually updates values of the stream function as they become available. The method 
is convenient to program and requires storage only for L number of points. However, 
this method was found to converge slowly and required many marches through the field 
to obtain an accurate solution. 

The method selected for this paper employs the Gauss -Seidel overrelaxation 
scheme given by 




R+1 


R 


Id 


= 


+ V 






R 


Id 


(1 ^ < 2) (6) 


Utilizing overrelaxation was found to reduce the number, of marches through the 
field by an order of magnitude over the Gauss-Seidel iteration scheme {v = 1), as indi- 
cated in reference 7 when applied to the Laplace equation. 

In the case of a profile where the surface stream function is known a priori (sym- 
metric airfoil and closed body of revolution), the previous equations are sufficient to 
determine the solution anywhere in the field. For the cases where the stream function 
ahead of and on the surface is unknown (cambered airfoil, airfoil in ground effect, or 
nacelle), an additional equation is required. The equation is obtained from the Kutta con- 
dition and it is assumed that the flow leaves the trailing edge alined at the average 
trailing-edge slope. For the trailing edge. 


- — = tan 6 
u 

By usir^ equations (2) 
d\jy/dr 


by usii^ the appropriate algorithms 




or by rearrai^ing 


'^l.r 3 i - J - fr •“ i'l'i.ni - \i-l) 


( 7 ) 
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Equation (7) applies only to the trailir^ e^e. Since the profile is a streamline 

itself, all profile points can be set equal to / 4 /. A and can be continually 

\ ’"trailing edge 

updated after each sweep through the field. 

The grid system incorporated in the present solution can have only one horizontal 
grid line intersecting the profile. (See fig. 1.) As a result of the simplified grid system 
employed, this present method is applicable only to slender profiles at an angle of attack 
of 0°. Naturally, one would not expect accurate results in the region of a blunt leading 
edge. 

The pressure distribution on various profiles is presented in the next section and 
is compared with classical solutions. The pressure coefficient can be calculated by the 
following relations: 

Cp = 1 - (u2 + v2) 


where 


u = 


J_M 

fK dr 



surface 


A listing of the computer program used to calculate the field stream function and 
surface pressures on slender profiles is presented in appendix B. Computational time 
for obtaining results with this program has been estimated at 1 minute per case for 
3000 field points and 200 field sweeps on the Control Data series 6600 computer system 
at the Langley Research Center. 

RESULTS AND DISCUSSION 


The present method has been used to compute the flow field and pressure distribu- 
tion on various profiles. The input field size was enlarged until surface pressures did 
not vary. Comparisons of results are made with classical or currently used techniques. 
All cases computed by the present method have required less than 1 minute computational 
time on the CDC 6600 computer. Computer storage is less than 50 000 octal for all 
cases. The detailed inputs are presented in appendix B for all cases which were run. 
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Airfoil With Thickness and Camber 

Figure 2 shows the computational grid for the NACA M6 airfoil which has 
12-percent-chord thickness and 5. 4-percent- chord camber. There are 2091 grid points 
in the field, the field extends 1 chord length upstream and downstream, and 2 chord 
lengths laterally. 

Figure 3 shows the pressure distribution on the upper and lower surfaces. The 
solid line represents the solution obtained by using exact potential theory, (See ref. 8.) 
The symbols show the results obtained by using the present method. In general, agree- 
ment is better on the upper surface where there are more grid points near the surface. 
(See fig. 2.) Agreement is poor near the blunt leading edge where the simplified grid 
system provides a coarse grid. 


Symmetric Airfoil in Wall Effect 

Figure 4 shows the computational grid for a symmetric airfoil whose geometry is 
given by the following relationships: 


r = 0.025 


r = 0.025 


1 - (1 - 2z)3 

1 - (2z - 1) 




1/3 


(0 ^0,5) 


(0.5 iz £ 1) 


The airfoil has a thickness of 5-percent chord. The wall is placed at 30-percent chord 
in the lateral direction. There are 2193 grid points in the field. 

Figure 5 shows the pressure distribution on the upper and lower surfaces. The 
solid lines represent the solution using a distribution of surface sir^ularities (ref. 9). 

The symbols represent the results obtained by usii^ the present method. The agreement 
is very good for this thinner airfoil except at the blunt leading edge where the grid spac- 
ing is too coarse to provide accurate results. 


Parabolic Body of Revolution 

Figure 6 shows the computational grid for a parabolic body of revolution with 
10-percent- chord thickness. The lateral field extends 1,2 chord lengths with 1071 grid 
points. 

Figure 7 shows the pressure distribution on the body surface. The solid line repre- 
sents the solution obtained by usii^ a distribution of surface singularities. (See ref. 9.) 
The symbols represent results obtained by using the present method. Once again, the 
agreement is very good in regions where there are grid points close to the surface. 
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Camber- Line Nacelle 


Figure 8 shows the computational grid for a camber- line nacelle which has 
2.135-percent-chord camber and a 13.671-percent-chord capture radius. The lateral 
field extends 9 capture radii laterally. 

Figure 9 shows the pressure distribution on the nacelle surface. The circular 
symbols are the results obtained by using the streamtube curvature method (ref. 4). 

The present method, given by the square symbols, shows very good agreement. 

NACA 1-89-100 Nacelle 

Figure 10 shows the computational grid and geometry of the NACA 1-89-100 nacelle. 
(See ref. 10.) Experimental data for this nacelle was obtained at a Mach number of 0.6. 

A translating center body was used in the tests to obtain results for different mass-flow 
ratios. 

The three-dimensional ’’Gothert" compressibility correction is incorporated in the 
present method. Because of the simplified grid system the center-body geometry could 
not be represented. However, results can still be obtained for different mass-flow ratios 
by deleting the Kutta condition. The Kutta condition or flow direction condition at the 
trailing edge sets the value of the surface stream function. If the mass flow is known, 
then the surface stream function Is known by 

^surface _ _m_ 

(See fig. 1.) Since the correct internal geometry is not used, the internal pressure dis- 
tribution will be invalid. However, the pressure distribution on the outer surface should 
be correct as long as the proper mass-flow ratio is used. 

Figure 11 shows the pressure distribution on the nacelle outer surface at a Mach 
number of 0.6 and a mass-flow ratio of 0.787. The circular symbols are the experimen- 
tal data. The square symbols are the values obtained by the present method. The agree- 
ment is very good except in the region of the trailing edge, where the proper flow direction 
because of the presence of the center body is not represented. The example illustrates 
the versatility of the present method to solve problems which are difficult to handle with 
classical methods. 
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CONCLUDING REMARKS 


A rapid numerical solution has been presented for the incompressible flow over 
thin planar and axisymmetric profiles at an angle of attack of 0^. The method uses a 
finite -difference field solution to the incompressible equation with a Gauss-Seidel succes- 
sive overrelaxation scheme. In spite of a simplified grid system and a small number of 
Surface points, results were in very good agreement with classical solutions for a variety 
of cases. The field-point formulation using the stream function is convenient for extend- 
ing the present capability to subsonic nonpotential inviscid and viscid flows. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 24, 1973. 
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APPENDIX A 


FINITE-DIFFERENCE ALGORITHMS FOR DERIVATIVES 
WITH UNEQUAL SPACING 

For second-order central difference approximations to derivatives of a function 
across an interval: 


i-1 ah i 


h i+1 


the function is of the form 


= 2Lq + SL^Z + 3.2^-^ 


because 


aj + 2a2Z 




Sketch (a) 


'^ i -1 ■ ^0 


(z = 0) 


i^i = + a^ah + a2a\^ 


(z = ah) 


'^i+1 = "(^i-i + (1 + a)haj^ + (1 + a)^h^a2 


(z = (1 + a)h) 


where ah and h are relative point spacings. (See sketch (a).) As a result, 


ah^(l 
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APPENDIX A - Concluded 


Therefore 


aj = 


1 

ah(l + a) 


+ (1 + - (1 + 2o!)i//. 




+ (1 + Q,)(l - Of),;, - 1 

dx Ij ah(l + a)L ^ ^ ^ iJ 


d^ 

dx^ 


ah {1 + a) 


“’^i+l ■ *^i-l 

+ a) L J 


H 



APPENDIX B 


COMPUTER PROGRAM FOR CALCULATING THE PRESSURE 
DISTRIBUTION AROUND SLENDER PROFILES 

The calculation procedure described in the main body of the paper for obtaining 
pressure distribution around slender profiles has been programed for high-speed digital 
computation. The purpose of this appendix is to provide a description of the necessary 
input and available output as well as a FORTRAN IV (ref. 11) listing of the source pro- 
gram. Input listings for the profiles studied in the main body of the paper are presented. 
The output listing for the cambered airfoil is presented as an example. 

Description of Program 

The program reads the profile coordinates and field geometry. The grid is gen- 
erated and initially labeled as uniform flow. The program marches through the field 
updating the value of stream function at each point accordir^ to the differential equation. 
At the trailing edge the Kutta condition is used instead of the differential equation. When 
the field sweep is completed, all profile points are given the same value of stream func- 
tion as the trailing edge, and a new sweep is started. The program proceeds for 
200 sweeps through the field, at which time the surface pressure coefficients are com- 
puted with the latest values of stream function. The results are then printed. Computa- 
tional time for the case of 3000 field points and 200 sweeps through the field is 1 min- 
ute. The required storage is 50 000 octal on the CDC 6600 computing system at the 
Langley Research Center. 


Program Listing 

The FORTRAN IV listing of the source program used on the Control Data series 
6600 computer system at the Langley Research Center is as follows: 


DIMENSION XOLOl) «Y0LT<3H fYOL0(3I» *Z < 1011 .9 (ftl » *PSI ( 1 01 » fil » » 
ISLOPETOn ♦SLOPE0<3H ,DELZ(101» .iJTAa(4n 
REAL MACH, mass 

NAMELIST/NUM/x 0L,Y0LT,Y0LR*R80L«N)(*SCALEZ.SCALER» JTE,r>ELZS»NC* an 
1 »MACHiMaS5 
500 READ(5»NUMI 

1F(E0F.5> 200,99 
99 WRITE<6,NUM) 

P 1=3. 141592653589793 
RB0L=RR0L*S0RT (1 .-MACH**2) 

00 1 1=1, NX 

Y0LT(II=Y0LT<n*SQRT<l.-MACH**2t 
1 YOLB ( I > =Y0LR < 1 1 *S0RT ( 1 .-M ACH**?> 

5CALEZ=SCALEZ..1E-0R 
SCALER=SCaLER*.1E-0R 
C GENERATION OF COMPUTATIONAL GRID 
AL=1. 
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APPENDIX B - Continued 


RB=RBOL*AL 

IF<JTE.EQ.1IOFLR=.060 
1F(JTE.F0.11G0 to 17 
DELR=«80L/( JTE-1 > 

17 CONTINUE 

II=?*NX-1*IFIX( <2.«SC4LEZ-l.)/nEL7S) 

ILE=(I1-NXI/2+1 

ITE=ILE*NX-1 

JJ=TFIX f ( 1 .♦scaler >*«R/DFLR t ♦! 

IFUTE.EO.l) JJ=IflXISCALF»/OFLR)^l 

R(n= 0 . 

DO ? 

? R[J)=R(J-1 I+DELR 
IMID= (NX-1 I /2+1 
DO 30 IP=?*IMI0 
I=TLE-1B*I 

OEL7(I)=XOL(lR)-XOL(lR-n 

30 CONTINUE 
IMF=ILF-(NX-U/2-I 
00 31 T=1 ,IMF 
0EL7(I)=0ELZS 

31 COMTINUF 
IMR=TTF+ (NX-1 )/2 
ILTM=IMB-1 

00 3? I=ITt«ILl8 

ib=ite-ile-(I-ite> 

OELZ (It =XOL (IB>l)-XnL(IR) 

3? CONTINUE 

00 33 l=lMfl.II 
0ELZ(I)=0EL7S 
33 COMTINUF 

k=itf-i 

00 3S I=ILC»X 

35 OELZ (It =XOL ( I-!LF+2i-X0L (l-ILP+l) 

Z(lt=0. 

00 34 1=2,11 

,34 Z(I)=Z(T-J)*DFLZn-l) 

C NUMBERlMfj OF POUNnARY STREAM FUNCTION 
00 3 1 = 1 . 1 I 
PST (I,?)=SOO. 

3 PST (l.n= 0 . 

DO 4 

PSI (1 ,.)) =PSI (I ,2)»<R(,M/DELR) **NC 
DO 4 I=?,Il 

4 PSI (I,J)=PSI(I-U J) 

C FIELD MARCHTN(3 LOOPS 

K=TI-1 
L=JJ-1 
IC(HJNT = 0 

10 IC0UNT = IC(3UNT*1 
DO 5 J=?»I. 

no fi 1 = 2 , K 
A=1 . 

IF ( ( J.FO. JTE ) .ANn. (T.OE.ILE t .aNO. (T.LT.ITF nooTn^v 
IF( (J.FO.JTEI .and. (I.FO.TTF) too TO 11 
GO TO 12 

11 IF(MASS.GT.10.IGO to 15 
PSI(I,J>=PSI{1 ,JTE)*MaSS 
GO TO h 

C KUTTA CONDITION CALCULATION 

15 C= ( YOLT (NX-1 ) -YOlP (NX-1 ) 1 / (2,»{ XOL (NX) -XOL (NX-m t 

B=DFL7 ( !♦! ) /OELZ ( 1 ) 

PSI (I , J> = ( (1 .♦81 * < 1 .♦«) "'’ST ( !♦ 1 , I) -PS I ( I + ?,.J) -C"P"nFLX ( I ) * ( 1 .♦5> 
1" (PSI ( I , J^l 1 -PSI ( [ ,.)-l ) ) /(?.*OFLP) ) /(RorP.^Pt ) 

IF( (J.FO. jTE) . and. (T.FO.ITE) too TO G 
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APPENDIX B - Continued 


1? IFr ( J*rO* ( JTF+D ) n*r,F*TLF) .ANO. 4T»Le*TTFMGn TO 

IF ( ( J.fO*(JTE-l> » .ANn* (T*r,F,ILE) •AND. U*LE,!Tfm 
1A=U-Y0LB(I-TLF^1 \ /OELR 

H = 0'^L7(T)/0FL7(I-I \ 

FSJ ( I , J> =PSJ ( T f *AN^ ( 

6 (PSn I * J-l U •♦A) <»nFLR<^nFLR) fNC-1 > / f M A> <>OFLP*P ( J1 > y ♦ 

IPST ( I • J+1 > ^ (2./ < { 1 •♦A)'»nELP»nELR) - (NC-l ) / f A»( i , + A > »nFLP*P < J ^ > i ♦ 

2 (P«P5I + *OFLZ ij~}) ^Hfl? ( t-)) U ) / 

7 {2*/(R«nEL^ ( I-l) <»nFLZ ( T-*l ) 1 ♦ ?• / ( a^^OFUR 
3*DELP)-(1*-A)^(Nc-I)/( A« nELR*P( ) ) 

no) > 

GO TO F 

R A = K-Y0LT(I-ILF*1 >/nFLR 

R=HFLZ (T) /OEL7(I-I > 

PST Uf J)=PSI (I^J) *AN*( 

6 (PST ( 1 1 J-1 » * (?• / ( A*nFLR*nFLR'» ( l , *A) > + (HC-1 > / ( w < j> *a»dflp 

1* ( I *+A> U ♦Psr ( I f <2,/ (OFLPOOFLP* n .♦A) > -A<> (MC-1 ) /{p ( J> <*DFLP 

5 < 1 , ♦ A > ) ) ♦ 

? (4*»PSI <T-TiJUPSTU^UJ)lo<2*/(P»M. *p ) h>OFL Z ( T - i > ^’■nFLZC I-l ) ) ) t / 
7(2./(rooel^( i-i> *on.zn-n ) +2, / (A<>nELp*nFLP)-< a-ui i 
3/ (P{ J) <>A»nFLR> > 

A-PST no) > 
ff. continue 
S CONTINUE 

00 9 I=ILF»ITE 
9 PSn I, JTF) =PSI UTE* JTF) 

IF (icOUNT,EO*^0i/) GO TO 20 
GO TO 10 
20 CONTINUE 
IT=NX-! 

DO 14 I=?*IT 

SLOPFT M ) = ( volt < !♦ 1 > -VOLTU-1 M / ^ XOL ( T ♦ I > - XOL I T -1 > i 
14 SLOPFRU )=-!•» <VOLR(I»1 )-yOLR (T-1 > ) / ( xOL U ♦ I ) -XOL ( T -1 ) ) 

SLOPFTd )=Y0LT<2)/X0L(?) 

SLOPER ( U=-1.*Y0LR (?) /XOL (?) 

SLOPE! (NX) =-l *«YOLT (Nx-P / ( XOL (NX) -KOL (NX- 1 > ) 

SLOPER (NX) =Y0LB (NX-1 > / ( XOL ( NX ) -XOL ( NX- 1 ) ) 

UOspSl U f3)/(?.<>P(2)«‘*<NC-l)*nFLP) 

It=ILE 

IT=ITE 

IF(JTF.E0.n IL=ILF*1 
IF(JTE*FO.n IT=ITf-l 

C calculation or velocity COMPONFNTS and PPESSURF COErFTCTENT 
00 13 I=ILiIT 

A-1 ,-YOLT ( I-ILF* I >/r>FLR 

UZT = l . / ( (P ( JTf) ♦ YOLT ( I-ILE + l > ) (MC-1 ) «A*OEL^^n • ♦ A) ) 

1* (-1 ,'»A**?*PSI ( 1 ♦ JTE^?) + ( 1 •♦A) »<>?«PSI ( T • JTE+n-( 1 
2PST <r,JTE) )/UO 

URT=UZT*SLOPET(I-lLE^l ) 

UZT= (UZT-l • ) / ( I •-MACH*«?) ♦ 1 * 

URT = URT/S0RT < 1 

CP T=1 •- (U7!**?*UrT**2) 

IF (MACH,GT.O#) CPT= ( ( 1 •♦WACH*»2*CPT/S. > *«3*5-l* ) / ( . 7 «^mach**2) 

IF(JTE.EO*1)GO to IR 

A-1 .-YOLB ( I-ILF+1 ) /HELP 

UZp=I , / ( (R( JTf) -YOLR ( T-ILE*D ) «* (NC-l) (1 •♦A) ) * 

1 ( (U + 2**A)*P5I (KJTe>“(U+A)<^*?«PSI (!♦ JTE-1) ♦A*<»2*PSI (1^ JTE-2) ) /UO 
URBsUZR^SLOPER (I-ILF^ l ) 

UZb= (UZP-l •) / (I .-MaCH^*?) -fl, 

URB=URB/SqRT < 1,-mach«^»2) 

CPB=1,-<UZB**?+UPR<m^2) 

IE (MACH.GT . •! ) CPR= ( ( I • ♦MACH*-»p*CPf' * ) / ( .7»MACH*«?) 

18 K=I-TLE*1 

WRITE (6, 400 > XOL (K) ♦CPTtCPR 

AOO format ( 2X4HX0L=El2#A#?X4HCPT-E12*4»2X4HCPRs:ei2*4) 

13 continue 
GO TO 500 
200 STOP 
ENO 
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APPENDIX B - Continued 


Input 

A single case consists of the pressure distribution around a given profile. It is 
necessary to input the profile and field geometry, planar or axisymmetric code number, 
value of over relaxation factor, Mach number, and mass-flow ratio. For the loading rou- 
tine used in the program, any column except the first may be used on the input cards. A 
decimal format is used for the input quantities. A description of the required inputs and 
fortran variables used by the source program are as follows: 

fortran ^ 

variable ' Description 


$NUM 

X0L 

Y0LT 

Y0LB 

RB0L 

NX 

SC ALE Z 
SCALER 
JTE 
DELZS 

NC 

AN 

MACH 


arbitrary name required by the loading routine to define the input data block 

array name for axial profile coordinates 

array name for upper surface coordinates 

array name for lower surface coordinates 

capture radius in percent of chord 

number of surface points input 

axial extent of field in percent of chord 

lateral extent of field in percent of capture radius 

jth value of approach grid line 

standard axial increment for far field. Multiple of largest surface axial 
increment. Must be compatible with SCALEZ 

code number for planar or axisymmetric 
NC = 1 planar 
NC = 2 axisymmetric 

over relaxation factor 
1 i AN < 2 

free -stream Mach number (use for axisymmetric only) 
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APPENDIX B — Continued 


FORTRAN „ . 

variable Description 

MASS mass -flow ratio 

$ denotes end of case (column 2) 


The system loading subroutine in the program (namelist) is very flexible in that the 
order of the input cards is unimportant and successive cases can be run by repeating the 
identification and $NUM cards followed by only the changed parameters and a $card. 
Experience has shown that 1.85 to 1.95 is a suitable value for the overrelaxation fac- 
tor AN. The input listing for each of the profiles is given. 


Input for NACA M-6 Airfoil 


XOL 


YnLT 


YOLB 


RrOl 

NX 

SCALBZ 

ScarCR' 

JTE 

delzs 

NC 

An 

MAy 
Mass' ‘ 
send 


= 0.0* O.IBSE-Ol, 0.?SE-(H* O.Sr-Ol, 0.7S^-01. 

O.ISE + OO. + 0.AE'*’0o. O.bE + 00* 

0.7E*0(1* 0.8E + 00* n.BSr + OO* (i.BE + Oo. 0,9?;SF*00. 

o.byse^oo. C.bayRe+oo* n.iE^oi. r. t. t. i. t* t* t. t. t* t. 

= 6.0* 04197E-01. O.RBlE-01. 0.403E-O1. d,494F-01, 0.S71F-0I. 

6Vf'fi?F-0i ’ 0.7">sF-nr. ‘ o!i;8??F-0'l . * O.BOSFr-01* r 

0.,603F-'01. 0.30FF-01, 0.B3F-01. O.iSBF-Ol, 

0.1??E-0l. 0»fl«F-0?» 0.44F-02* 0.??F-i1?t 0.0. T. T* T* T« T* I. 

U T* I* I. 

= 0.0* 0.17(^E-0l. 0.??E-01* f).?71F-0l. Oi1Cl3P-Ol. r.3?4F-01. 

0.347F-01. 0»36?E-0l» 0.37OF-01. 0.39F-01. 0.3B4F-01. 

■‘o.3B2e-oi, o.34BF-oi. ■ o.?s3F-o"r.”' or?3F-"6i“ ir.'TTyF^r; 

0.143E-Oli a.lOBE-01. O.S9F-0?. 0 ,?QE-02 . ' 0 .0 . T . I . Tt I. 1 . "l . 

It T* I* T. 

= O.BF+Ol* 

= 21 . 

= O.lFtOl* 

■=■■ oVirtoii ■ ' “ - ' 

= 21. 

= 0. IFtOO* 

= 1 * 

= O.lBSEtOl. 

= 0.0. 

= 'b.BEi'o?* ‘ ■" “ 


Input for Airfoil in Wall Effect 


$MUM 

XOL 


YOLT 


yolr’^ 


0,0f 0.1?SF-UU 0,S^-*0U ..0*7SE-Olf O^IF + OO* 

0.15F + 00* + o.4F + ri()4 0'i5F+b0t_ .o/f^F + OO* 

0.7F + 00* 0*8E + 00» 0*QSF + 00% 0*RF + 00» n.QRF + OO* ~ 

0*97FE^0O? O.lc+Olv U U It' It T. T, ‘Vt T, 

0# \ 04SSQ R? fl}S0l4f^n ] » j 30 61^6^1 1 241SF-01 1 

0-1617H1R4067401F-‘01 . 0 * 1 8P?3 07R145 1 704F-0 T t 0 • 1 ^568 ?4aSQ 1 S5 1 l_F-0 1 t 

0.3173143Q633S9PF-0] t 0 * ?3QS? 1 B1 I . 0., ?44F4866?3?SA4^0 1 t 

? * f ^ ,._pt P?5F-0 1 V 0.*?34£-01^ O.IP/SF-Ol* 

0 • 1 Fa^?5E-01 1 0,l2PE-01t O ♦ Q646R74;4‘^QQOOQE-"02t. 

0 . 677500000000 r>?F-0? • 0 . 3‘^^S6?4^9^9'?Q9F-Q2_t 0 ■ 1 8PRR 1 , 

. ^±S) Xi_l t It It I t 1^ L._ Lt___ ' 

OtOt 0. iri4'5505?RiF034E-bl^ 0 . 1 306 1 1 24 1 5F-0 1 * ' _ 

Ot l6178lP406740IF-*0,It 0 , 1 R?007334S1 704E-0 1 t .0 * J Qft8?4RS9.1 5Sj lF-0 1 

0-?173343Q613S9RF"'01 • 0 • 230S? 1 R 1 460 P9PF^0j t ' 244F4B66?.32S44F-0 1 . 

0*?4933lS476l 193F-^01 t 0*P5F-01t 0.?4RF-01.t 0*234E-0lt <^*l96F-0lt 

0. 1 642yr-0 I ’ 0#122E-ult 0 •96468749QgqQg9E-02t 
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APPENDIX B - Continued 


Q 7750 o.o^Q.Q.s>J o?r-*o?^ 0 * 

0*04 T* It T» !♦ ?.♦ I* T* T* !•_ 


. rbol 

= 

0.3E+00i 

NX 

:: 

21 4 

scaleh 

= 

0. ] OOOOOOOOlE+01 4 

scaler 

= 

0.^)000000(I01F*01. 

JTE 

= 

74 

DEL7S 


0.1E*00* 

NC 

= 

l4 

AN 

= 

0.19FF+01 4 

MACH 

=: 

0.04 

mass 

$FMD 


0.?E4024 


o;ooo 

o.oob 

0.000 

.010 

.010 

.on 

.013 

.013 

.026 

.016 

.016 

.050 

*018 

.018 

.075 

.0^0 

• o?o 

•100^ 

• 0?2 

.0?^ 

.isn 

.0?3 

.023 

.200 

.0?^* 

.0?^ 

»J0 0 

.o?s 

.0?F 

.^00 

.02s 

• 025 

.500 

;02S 

4 025. 

.600 

.023 

.023 

.700 

.o?o 

.020 

.BOO 

.016 

.016 

.850 

.012 

.01? 

.900 

.010 

• 016 

.925 

.0Q7. 

,...•0.07 

.960 

.0 0^ 

.6oi 

.975 

.002 

.00? 

.9br 

o;ooo 

0.000 

1.000 


o~ i R?Q S 1 S6? 4qQ90_r_"n ? 1 . 





Input for Parabolic Body of Revolution 


$NUM 

XOL 


YOLT 


yoLr 


0 

0 * 

0 . 

0 . 

_0 

0 

0 . 

0 * 

0 . 

I 4 

0 

"" 0 . 

0 * 

0 . 


04 0 * 1 ^>SF- 01 * 0*?5F-Olt O.RF- 01 *_ 

0 .?F + 00 t 0.3F + 00 * 0*4F + 00 4_ + 

7E4-004 O.BF-400* f?,RSF' + 00 * _ 0»Qr.4*0 0 _* 

91SF.*^6* U*9B75F + 004 0 • 1 1 4 T 4 J ^ T • J^4_J 4 T 4 T 4_J V_[ v. T 4 ^^ ^ — 

0 * Q*?^ 6 B 75 F- 0?4 6*4R7SF-0 ?4 0 *^SE- 0 ?< 0 . V3R75^-01_?_ 0 • 1 RF- 01 a_ 


0 ,SF- 0 l 4 


pl5^0W~ 0*3:?F-05 4 0*4?F-0l4 j_4 

4RF-01 4 0 4 0,3] QQQqQQQQ QQQ Qt^-g 1 4 

0 . 1 ? 0 *‘^44^Q990Q9Qqg6F-n?4 

4874<99qOq999Q7F“-0?4 0 • ?46B749999Q9Q 1 F -02 4 .0*04 T 4 1 

p! ^0*176 Fj^ 1^ 0_^22F^l *_p.^ 0*363F-0l4 ri,3?4E-01 

3 4 7 F - b 1 * ^ 36 ?F - C 1 4 0*37 9f -0 1 4 

382 F-014 0*348E^Ol4 0.?R3F-014 

l43E-0lt 0,108E-014 0,Fqr-0p4 fi.P9F-ft24 O-O. t* 1% I 4 ! 

T4 I* T4 


0,?SSF-0 1 4 


it ?4 1 4 T . 


0 * 39 F- 0 l 4 

0 . 23 F- 6 i 4 

a,? 9 F -024 


■ 0 Vl 77 F-bl 4 

0.0 4 T 4 1 4 1 4 ! 4 i 4 T 4 


rbol 

NX 

— 

0.01 

21 . 

scalez 


b. lOoooooboiE+oi 

scaler 


0. lPOOOOOOOlE+01 

jte 

= 

l4 

oelzs 

- 

0.1F400* 

NC 

= 

2 4 

AN 

= 

0.195F + 01 4 

MACH 


0.04 

mass 

JFND 

~ 

0.04 
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APPENDIX B — Continued 


0.000 

0.000 

0.000 





♦ 00? 

.01 A 

.013 





♦ CDS 

.0?? 

.n?B 





.009 

• 6^7 

.050 





.01^ 

i030 

.07B 





.OlR 

.03? 

*i00 





♦ 0?5 

i03S 

.ISO 





.03? 

.036 

.200 




— 

;04? 

;03A 

.300 





:048 

• 03Q 

.400 





;050 

.03Q 

.500 




~ • 

.04P 


.600 




— 

.04? 

.035 

.700 




■ 

.03? 

♦ 0?R 

.ROO^ 






;o?5 

„ _s.0?3_ 

jB50 





;oi8 

.018 

.900 




- „ 

;oi4 

.014 

,9?5 




■ ' ' — 

.009 

.011 

.950 






.005 

.006 

.975 




‘ ■ ' — 

;oo? 

.003 

.9RR 




— - 

OiOOO 

JLiOM. 

1.000 




- • ’ 


36 SI ?i 

1 iiooo 
? iiooo 


Input for Camber- line Nacelle 


snum 

XOL = 0.0* 0.1P5F-01* 0.?SF-01« O.S*^-OU 0.75F-01* O.lF + 00. 

0.15F*00. 0.?F*00. C.3F*0Oi O.^F+OO* O.SF+fiO. O.^E+ho* 

0.7F*00i 0.3E«00* 0.ftSF»00. 0.ciF»OO, 0.'^?SF-K10» 0.'3SF*0O, 

0.<77FF»00. 0,9075F»00» 0.1F*0l. !. T. 1. T» I. T» !• I. !• T, 

'N?7F-o?. o.^l^Fj^;^ o.in7aK-ni , 6. i ?<=; sf-oi . ' o.i407f-oi. 

0.16£tF-0i, 0.1794E-6r« 0.i99'4F-01%’ 0.>101F-0!, 0 . ? nSP'-O f. 

0.210IF-01* C.1994F-01* 0.1794F-01. 0.164F-01, 0.1407F-01. 

O.iasSF-OJ* 0,107 «f- 01, 0.913F-0?. 0.37e-0?. O.o* T* 1* r« !• I" 

!• I* I* T, I* 

YOLB = 0.0* -0.B7F-0?* -0.B13F-0?* -0.107BF-fH, -0 . ) ?<sSF-0 1 . -O.UnrP-Ol. 

. . -Of 1794F_-01 , -0.1994F^1_^ -q^iOlF-OJ_. -0_.213BF-iJJ . 

-0.2101F-‘0!i -0.1994E-01 • -0 . 1 794E-0 J ^ -0. Ife^F^o i " , j 4 q 7I- lo j , 

-0.125BE-01* -0. 1 07BF-D1 , -O.413F-0F. -0.37e-C2-. O.O* I* !• t. I. T* 
T* I* li !t It 
RbOL = 0.13671F*00. 

NX = 21 . 

ScALEZ ^ O.IF+Ol* 

scaler = 0.9E + OU -- --- - - 

JTF = B, 

DELZS = O.BF-014 

NC =2* 

AN = 0.1BBF*01. 

MaCH = 0.0* 

MASS = O.PE + 0?* 

SFNB 


Input for NACA 1-89-100 Nacelle 


$NUM 

XOL' 


iQki 


L Q.i2^E-(fi. nc?sE-oT« *T.^-irr* OTTE^oiri ffvr^vw* 

!_ 0.3g *60 * 0.3sF*oo, ~0.4r* 00* 0.4SE*g07-^'.3FFOoT 

® • -."zf rP..P_f_. P . 6E +0 (T.'6^ oo t 67Wi 0 0 » (TTTBr^Wi — tr. ftE*0o '.~ — 

P ^9SE*00, 0.9E.OO, 0.95E-*od, 0.975F*00» 07^975 riW7^ir.TF*6T~. 

■ - .P---9.f.-_,i;.Ai3E7_p 2, 0.604E-62* 0.863E-02* 0.1229E-01, 67T4BfE^oTi 

__P_..1 908«=--01* O.l fteSE-0 1 . 0 . 1904E-01 ♦ " oVm^E^*. 

— P f i ?P.^F --'Ll!. _P . 1 91AF^ 1. 0.I904E-01* 0.f904E>01 * O'* T^^F-TT ♦ — 
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APPENDIX B - Concluded 


0,L<30J^-Ql» 0 .iBQ4r-0lf 0.179 6E-01 . 0.1577li-01f 0 f^l227E-U.! 

0*717E-0P» 0-179E-^* 

YOLr = • - 0 ;i 9 E- 03 r- 0 . 4 F-b 3 \ -Q. 65 F*Q 3 ^ * 0 , 16 ~ 3 F- Qg ^ ^Q^gTr^ Ogf 

_0,/V46>*Q?* -0*63fiE-02f -0,75E*02# -0~.7 5E-D?» **0«75E*Qg» 

-oi 7 s_r-o’?^ -o# 75 Fj^^^_i^o* 75 EjiLlilr®i!^ 

_q . 7'^F-J? ♦ ~ *0V74ftE- 0^? ♦ «b>^lF- 02< --Q* I52E->02 • - 6 


RBdl' 

NX 

= 

0*1SA04E4-06» 

2S, 

scftCff 

= 

O.JP*01t 

SCAt ER 


0.9F*01f 

JTE 

= 

Pf 

delzs' 

= ■■ 

O.lF+00, 

NC 


?♦ 

AN '"■ ■ 

= “' 

''d.iA5E*01» 

MACH ■“ 

~” 

0.6F^00. .. 

M'ASS 


6.7B7E*00. 

$ENr> 




Output 


The output listing consists of the upper and lower surface pressure coefficient, A 
sample output listing is presented for the cambered airfoil case. 


XOL = 
X0L = 
X 0 L = 
XOL.^. 
X 0 L = 
X0L = 
XOLs 
XOL = 
XOL = 
XOL^. 
XOL = 
XOL = 
XOL = 
XOL = 
X0L = 
XOL = 
XOL = 
XOL = 
XOL = 
X0L = 
XOL = 


0. 

rPT= 

1 .?S00F-02 

CPT= 

2.5000E-0? 

CPT« 

5.0000F-0? 

rPT- 

7.5000F-02 

CPT^ 

1 .OOOOF-01 

rPT= 

1 .5000E-01 

CPT~ 

2.0000F-01 

CPT = 

3.000CE-01 

CPT = 

^.2000E-01, 

CRT- 

5.0OOOE-O1 

CPT = 

f>.OOOOE-Ol 

CPT = 

7.0000E-01 

CPT = 

e.ooooE-^oi 

CPT=t 

8.5000F-01 

CPT = 

Q.OOOOE-Ol 

.,CRT= 

9. 250 OF -01 

CPT= 

'i.SOOOF-Ol 

cpt= 

<5.7500F-01 

CPT= 

9.8750F-01 

CPT = 

l.OOOOF+00 



A67ftF.-fll 


-5.221 7F-01 

C^R = 

-6.5589F-0? 

rpR= 

-2.2225E-01 

tPM = 

-3.1975F-01 

CPR = 

-4.0fr71E-0l 

r^R- 

-5.???7E-nl 

CPR = 

-5.4801F-01 

CRR* 

-5,«;010E-0l 

CPR = 

-4.F875E-01 

C.°R= 

-3.5354F-6i 

CPP = 

-2.1705E-01 

CPR = 

-1 ,36IFF-f|l 

CPo = 

-7.n207E-03 

CPP = 

4.17805-0? 

CRR = 

9.4420F-02 

CPR = 

1 . i 3A>F-0i 

CPR = 

1 .3221E-01 

CPP=? 

1 .8229E-01 

C^R = 

2.H51F-01 

CPR = 

2.41Q8E-01 

CPR = 


-Q*i*9SlF-/ I 

-?.7.^4,0F-f?l 

- 3 . 1979 E- 0 ]. 

-3.riRP>^F-6i 

“?,90S4F-'»1 

-p.3Qa?F-f^i 

-1 •RSI IF-CJl 
- 1 *S 9 qRE -01 
-K 7272 F -01 
-UR 1 S 7 E- 0 I 
-KB0^4F-01 
-1 •s;^ 4 SF-rii 

-R* 467 FE-i')_? 

-3*i0^3F-0? 
-1 .PH^^F-CO 
R.3SR1F-02 
KOHPSF-01 
1 •S993F-01 
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Figure 4.- Computational grid for airfoil in wall effect. 
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Figure 6.- Computational grid for parabolic 



of revolution. 





Figure 8,- Computational grid for camber-line nacelle. 
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Figure 10 *- Computational grid for NACA 1 - 89-100 nacelle* 






